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Abstract 



We discuss in detail the physics reach of an experimental set-up where electron neutrinos (anti- 
neutrinos) produced in a beta-beam facility at Fermi National Accelerator Laboratory (FNAL) 
are sent, over a distance of L ~ 1300 km, to the Deep Underground Science and Engineering 
Laboratory (DUSEL). A 300 kt Water Cherenkov (WC) detector and a 50 kt Liquid Argon Time 
Projection Chamber (LArTPC) are considered as possible detector choices. We propose to use 
^'^Ne and '^He as source ions for i/g and fe beams respectively. The maximum Lorentz boost factor, 
7, available for these ions using the Tevatron are 7Ne = 585 and 7He = 350. This particular set-up 
provides the opportunity to probe the first oscillation maximum using beam and the second 
oscillation maximum using the i>e beam which helps to evade some parameter degeneracies. The 
resulting physics sensitivities for ^13, CP violation and the mass hierarchy are compared to those 
of a conventional superbeam from FNAL to DUSEL. 



*Electronic address: 3anjib.at.vt.edu 
^Electronic address: pahuber.at.vt.edu 



1. INTRODUCTION 



The discovery of neutrino mass is one of the first clear cases of physics beyond the Stan- 
dard Model. Neutrino oscillation is now firmly established as the leading mechanism for 
flavor changes in neutrinos This has triggered substantial interest in precision mea- 

surements of neutrino properties, with a special emphasis on oscillation experiments. A new 
generation of long baseline [§, |^ and reactor experiments is currently trying to measure 
6'i3. This is considered to be the first step to a large scale experimental program of long 
baseline neutrino experiments which aim at determining the neutrino mass ordering and to 
study leptonic CP violation, see e.g. ||Tl], O]. Due to experimental constraints on the avail- 



able neutrino flavors in both source and detector, the most promising oscillation channels 
are the ones where a muon (anti-)neutrino oscillates into an electron (anti-)neutrino or vice 
versa. In spite of using both neutrinos and anti-neutrinos, a serious problem with all long 
baseline experiments involving these channels, arises from discrete degeneracies which mani- 
fest themselves in three forms: the (^i3,5cp) intrinsic degeneracy |TB[, the {sgn{Amli) , 6cp) 

?23,7r/2-^23) 



degeneracy |I4 

m 



and the {623, n/2 — 623) degeneracy |T^. This leads to an eight-fold degen- 
eracy with several degenerate solutions in addition to the true one. The presence of 
these degenerate solutions can severely limit the sensitivity of an experiment. 

a flavor pure 



In this paper, we will consider a beta-beam, as originally proposed in [17 



beam of electron (anti-) neutrinos from the beta decay of ions of short lived isotopes, which 



are accelerated to a Lorentz factor, 7 > 100. 



and there is a vast literature on the many different options for such a facility [18-38 



This idea has attracted a large community 

We 



present a study of beta-beam facility at FNAL. A somewhat similar idea has been explored 
in [^: there, however the emphasis was on possible synergies between the new facility and 
NOt/A. We, on the other hand, focus on a full fledged, stand-alone beta-beam. We will use 
^^Ne (for i/g) and ^He (for z/g) ions. The maximum Lorentz boost factors using the Tevatron 
are 7Ne = 585 and 7He = 350 which yield energy spectra which peak at 2.3 GeV and 1.4 GeV, 
respectively. Our detector is located at DUSEL p9|-p[1 at a distance of L ~ 1300 km from 
FNAL. We propose to use a 300 kt WC detector or a 50 kt LArTPC as a possible detector 
candidate. The first and second oscillation maxima for the FNAL - DUSEL baseline are at 
2.5 GeV and 0.8 GeV for Amg^ = 2.4 • 10~^eV^. Therefore, this set-up provides an unique 
opportunity to work at the first oscillation maximum using the u^. beam and the peak energy 
of the z/e beam is very close to the second oscillation maximum. 

The paper is organized as follows. We begin with a brief description of the FNAL based 
beta-beam facility in Section 2. In Section 3, we deal with the relevant oscillation probabil- 
ities. In the following section (Section 4) we describe the characteristics of the WC detector 
and the LArTPC; also, we introduce the superbeam which will used for comparison. We 
also present the expected event rates for these two detectors. The details of our numerical 
technique and analysis procedure are presented in Section 5. In Section 6, we present our 
results and provide a summary. 



FERMILAB BASED BETA-BEAM 



The concept of a beta-beam was proposed in [|T7|: a pure, intense, coUimated beam of 



electron neutrinos or their antiparticles from the beta-decay of completely ionized unstable 
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Energy [GeV] 

FIG. 1: The un-oscillated beta-beam flux spectrum arriving at a detector placed at DUSEL at a distance 
of 1300 km from FNAL. The red solid line shows the v,, spectrum generated from ^^Ne with 7 — 585. The 
black dashed curve depicts the Ve spectrum originated from ^He with 7 — 350. The blue dot-dashed and 
the brown double-dot-dashed vertical lines show the locations of first and second oscillation maximum for 
the FNAL -DUSEL baseline corresponding to Am|i = 2.4 • lO^^eV^. 



radioactive ions circulating in a decay ring. The first stage in a beta-beam are protons 
with an energy of a few GeV which impinge a neutron spallation target. These neutrons, 
then interact in a secondary target to produce the desired unstable isotopes, in our case 
^^Ne and ^He (see, table |). Both are noble gases and, therefore can easily diffuse out of 
the secondary target, where they are collected, ionized and bunched. Now, they can be 
accelerated and be put into a storage ring with long straight sections. The decay of the 
highly boosted ions in the straight sections then yields an intense, well coUimated, flavor 
pure electron (anti-)neutrino beam of known flux and spectrum, see e.g. Feasibility 



of this proposal and its physics reach is being studied in great detail [|T£|. The set-up we 
propose will take full advantage of existing accelerator facilities at FNAL |3^. The main 
limitation is the maximum energy per nucleon the Tevatron can achieve. The maximum 7 
available for the ions considered here, using the Tevatron as accelerator, are 7Ne = 585 and 
7He = 350. We have also studied the possibility to use and '^Li. These two isotopes have 
much larger endpoint energies, of more than 10 MeV and thus would allow to have a peak 
energy for both neutrino and anti-neutrinos around the 1st oscillation maximum. However, 
since they would require a much lower 7, they would yield a considerably smaller neutrino 
flux. We have performed a full numerical study for these ions for the full range of available 
7 and found that their performance is always worse than for ^He and ^^Ne . In a similar 
fashion, we have optimized the 7 for ^He and ^^Ne and the values chosen above represent 
the optimum in terms of physics sensitivities at the given baseline. Since, ion production 
requires a proton power of only a few 50 kW, the currently available proton intensities may 
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Ion 


tl/2 (s) 


Eo (MeV) 


Total useful decays 


7 (Tevatron) 


Beam 


(GeV) 




1.67 


3.92 


5- (1.1- IQiS) 


585 




2.3 




0.81 


4.02 


5 • (2.9 • lO^S) 


350 




1.4 



TABLE I: 



The beta-decay parameters, half-life tx/2 and electron total end-point energy Eq are 



shown in the first two columns |44|. In the third column, we list the total number of useful ion 
decays considered in this work. The maximum 7 available for these ions using the Tevatron are 
mentioned in column four. The peak energies of the and spectrum in the lab frame are shown 
in the last column. 



be sufficient. A new decay ring of about 14.7km in circumference (assuming 5T magnetic 
field and 36% useful decay fraction), however, needs to be built to store the ions. 

While the shape of the beam spectrum depends on the end-point energy Eq and the 
Lorentz boost 7 of the parent ions, the flux normalization is controlled by the number of 
useful ion decays per year in one of the straight sections of the storage ring N^. Table | 
depicts the relevant details of the properties of these ions and their considered luminosities 
and the choices of 7. We have assumed 1.1 ■ 10^* [y^ and 2.9 ■ 10^^ (z/g) useful ion decays 
per year for ^®Ne and ^He ions respectively p5 |. 

Figure |1| shows the un-oscillated beta-beam flux reaching at a detector placed at the 
DUSEL site at a distance of 1300 km from FNAL. We see from this figure that the (z/g) 
spectrum peaks at 2.3(1.4) GeV. The first and second oscillation maxima for the FNAL - 
DUSEL baseline are at 2.5 GeV and 0.8 GeV for Am^^ = 2.4 ■ lO^^eV^ Therefore, the 
spectrum is well suited for the first oscillation maximum whereas the flux is sensitive to 
the second oscillation maximum. 



3. THE Pe^ OSCILLATION CHANNEL 

The simulation work presented in this paper is based on the full three flavor neutrino 
oscillation probabilities in matter, using the preliminary reference Earth model for the Earth 
matter density [^. However, to explain the nature of neutrino oscillations as a function 
of baseline and/or neutrino energy, it is crucial to use approximate analytic expression for 
Pe/i in matter p7|^9l, keeping terms only up to second order in the small quantities 6^13 and 
a = Ami/ Ami B 13 

9 r. 9 n sm'^\(l-A)A] ^ „ sin^(iA) 
Peu ^ sin^ 023 sin^ 2^13 ^ + cos^ ^23 sin^ 2^12 ^ 



To 



. . ^ sin A A sin 1 - A A 
± a sin 26^13 sin 26*12 sin 26*23 sin( A) — — - — sin 6cp 
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Detector Characteristics 


WC 

(Both fi^ k e^) 
(Only QE Sample) 


LArTPC 

(Both /i^ & e^) 
(QE & IE Sample) 


H 1 /H 1 1 /^T O 1 \ /I O C C 

riClUClcll iVidbb 


OUU KL 




-Cjiieigy J. Ill eoiioiu. 


WP- A /WP- R 


u.z ^^e V 


Detection Efficiency (6) 


OU/0 


oU/o 


Energy Resolution {5E) (GeV) 


0.085+0.05VE/GeV 


V.VoO-^V.yJD^y rjf lOI v^H/ odnipie 

0.085+0.2v/E/GeV for IE Sample 


Bin Size 


0.2 GeV 


0.2 GeV 


Background Rejection 


WC: A/WC: B 


10-3/10-^ 


Signal error (syst.) 


2.5% 


2.5% 


Background error (syst.) 


5% 


5% 



TABLE II: Detector characteristics used in the simulations. The bin size is kept fixed, while 
the number of bins is varied according to the maximum energy. We use two different simulation 
methods (WC: A & WC: B) to treat the backgrounds in WC detector. Details can be found in 
Section 4. 



/ . N sin(AA) sin (1 - A) A 

+ a sin 2^13 sin 2^12 sin 2^23 cos A — — ^ — ^ ^^cos^cp, 1 

A {I- A) 



where 

Amli L ^ A 

^31 



A = — 1^, A^—^, A = ±2V2GpN,E. (2) 
AE Ami 



Here, A is the matter potential, expressed in terms of the electron density Ne and the 
(anti)neutrino energy E; '+' sign refers to neutrinos whereas '— ' to anti-neutrinos. 



EVENT RATES AT THE DETECTORS AT DUSEL 



Currently, two technologies for large underground neutrino detection are considered for 
DUSEL: either a 300 kt WC detector or a 50 kt LArTPC 0. In the following, we 



will describe their different properties as far as needed for this work and we summarize 
them in table We do not consider backgrounds due to atmospheric neutrinos for either 
detector. The timing information of the ion bunches turns out to be sufficient to reduce 
these backgrounds down to an insignificant level, see, e.g. the appendix of |26| . 



A. Water Cherenkov Detector 

The water Cherenkov technology is well understood on a large scale and has demon- 
strated its excellent capability to distinguish muons from electrons. In a beta-beam, the 
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appearance signal are muons from charged current z/^ interactions; this has the advantage 
that they are easier to distinguish from neutral current (NC) events than electrons. Nonethe- 
less, NC events, in particular those involving one or several neutral pions, are problematic 
and they are the major source of background. Above the pion production threshold, the 
background level depends on how well neutral pions can be identified and distinguished from 
muons. Here, we consider only quasi-elastic (QE) charged current events for the appearance 
signal. We assume 80% detection efficiency, e, for both muon and electron QE events (see 
table H). More precisely, we should consider single ring events, i.e. for signal events, these 
events are characterized by only one charged particle being above Cherenkov threshold. For 
the NC background, this signature can be mimicked if the two photons from a 7r° decay 
are so close that the subsequent Cherenkov rings cannot be separated. This is more likely 
to happen for energetic tt*^ since the opening angle between the two photons is determined 
by the Lorentz 7 of the parent particle. Note, that identifying the quasi-elastic events with 
the single ring events is a reasonable approximation for the signal events. For the spectral 
analysis, we include Fermi-motion by term of a constant width of 85 MeV in the resolu- 



tion function. The energy resolution for the muon and electron is 5% of yE/GeV and the 
resulting width of the GauBian energy resolution function is the sum of both terms. 

As mentioned previously, for NC background events there is no simple, physically accurate 
approximation to their rate. Therefore, we will resort to two different phenomenological 
parametrization. The first one, method A, is based on an actual Super-K based Monte 
Carlo simulation, whereas the second one, method B, assumes an energy independent NC 
rejection efficiency. 

In method A, labeled as 'WC: A', we follow the results presented in where the authors 
study the performance of a water Cherenkov detector in a beta-beam using the current 
simulation and analysis tools developed for the Super-Kamiokande experiment. Results are 
derived for a baseline of 700 (130) km and taking 7 = 350 (100) for both ^^Ne and ^He ions. 
As far as the background is concerned, the main outcome of this study is, that the major 
background events in the search for appearance signal are indeed NC interactions involving 
pions. In a NC interaction, the outgoing neutrino carries a large and generally unknown 
fraction of the incoming neutrino energy. Therefore, those NC events which pass the single 
ring selection criteria tend to be reconstructed with an energy much lower than the true, 
incoming neutrino energy. It turns out, that, because of this, it is possible to maximize 
the signal to background ratio by imposing a cut in reconstructed energy. Unfortunately, 
in reference only results for 7 = 100 and 7 = 350 are presented, therefore we will 
extrapolate these results to the values of 7 of interest for this work. We assume that the 
average energy of a mis-reconstructed NC event, i.e. the ones which pass the muon single ring 
selection, is proportional to the incoming, true neutrino energy and that the proportionality 
constant does not change appreciable over the energy range considered. This intuition is 
inspired by the form of the differential NC cross section da/dy, where y = E^/E, with 
E being the energy of the incoming neutrino and Eg the energy transfered to the target. 
Therefore, we will assume that the energy cut which effectively removes the NC background 
is proportional to the average true neutrino event energy (E) nc 





j cP{E)Ear^c{E)dE 
! <P{E)aNc{E)dE ' 



(3) 
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where 4>{E) is the (anti-)neutrino beta-beam flux produced at the source and aNc{E) is the 
NC cross-section. We define a so called threshold factor, Tf which is different for neutrinos 
and anti-neutrinos, but independent of 7 

ET = Tf{E)Nc, (4) 

with Et being the threshold in reconstructed energy above which there is no NC background 
left. It is straightforward to compute {E)nc = 1.57 (1.69) GeV for ^^Ne (^He ) ions with 



7 = 350. From reference p2|, we find that Et = 1 GeV is sufficient for 7 = 350 to eliminate 



all the backgrounds for both neutrinos and anti-neutrinos. Using equation ^, we obtain 

Tf = 0.64 (0.59) for ^^Ne (^He ) . (5) 



Reference provides also results for 7 = 100 and we can use these to test our as- 
sumption that Tf is independent of 7. For 7 = 100, we have {E)j^c = 0.48 (0.5) GeV for 
^^Ne (^He ) ions. Now using the values of Tf given by equation ^, we get Et = 0.3 GeV for 
both ^^Ne & ^He ions which matches exactly with the value of Et obtained in for 7 = 100. 
This nicely demonstrates that the parametrization in terms of Tf and (E) tvc works over a 
reasonably large range of energies. With 7 = 585 for ^'^Ne ions we have (E) = 2.62 GeV. 
Now using the value Tf = 0.64 given in equation |^, we obtain Et = 1.7 GeV. Summarizing, 
for method A, we use a threshold energy of 1.7 (l)GeV for 7 = 585 (350) for both neutrinos 
and anti-neutrinos and assume that there is no NC background above Et- 

In the second method, labeled as 'WC: B', we take a threshold of 0.2 GeV, which is close 
the production threshold for muons and take the background to be a 10~^ of the NC current 
rate. The shape of this background is identical to (j){E)at^c{E) and for simplicity, we use 
the same energy resolution function to smear the NC backgrounds that we use for QE signal 
events. 

The impact of these two different simulation methods (WC: A & WC: B) is quite different 
while we calculate the sensitivity and we will discuss it in detail in Section 6. 

For our superbeam result which are shown for comparison only, we use the identical setup 
as in Section 10 of ET . 



B. Liquid Argon Time Projection Chamber 

A LArTPC works due to the fact, that by applying an electric field, free electrons created 
by the passage of an ionizing particle can be drifted over large distances, 0{m), without 
distortion. This allows, to obtain three projections of the particle track by just reading out 
the surface of the volume. For this reason, it seems feasible to build very large LArTPCs 
at a reasonable cost. The three projections of the track can be used to reconstruct the 
three dimensional path of the particle with an accuracy of a few mm. The currently largest 



LArTPC ever built, is the ICARUS T600 module with a mass of 600 1. It is recognized 
that scaling T600 by at least two orders of magnitude requires considerable R&D. 

Since this technology is still in its R&D phase, much less knowledge regarding its perfor- 



mance exists. Thus, what we present here, are essentially educated guesses [Q. We divide 



the signal event into samples of QE and inelastic (IE) events because energy reconstruction 
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works quite differently for these two event classes^. For QE events tliere are only very few 
tracks, typically the muon and the proton, thus it is feasible to perform a full kinematic 
analysis. Therefore, QE events will have a rather good energy resolution. IE events, on the 
other hand, will produce a large number of tracks, which the detector will not be able to 
separate and reconstruct individually. Most likely, only the muon track and a shower-like 
object can be identified, while the muon still is reconstructed quite well, only summary 
information will be available for the shower, e.g. total charge. Thus, the shower energy 
resolution will be much worse than the muon resolution. Therefore, we will use two different 



energy resolution functions for QE and IE events |]5^ 



SEqe{E) = (0.085 + 0.05 v'^/GeV) GeV for QE , (6) 



SEie{E) = 0.085 + 0.2 JE/GeV GeV for IE . (7) 



We use 80% detection efficiency, e, for both muon and electron QE and IE events (see table 
U). We calculate the sensitivity using two different values, 10~^ and 10""^, for the background 
rejection factor. The present status of the simulation study of the LArTPC at the DUSEL 



site can be found in Section 10 of 

As we divide the signal events into two parts, we divide the NC backgrounds into two 
parts as well. At a given energy E, the NC backgrounds, which are relevant to estimate the 
sensitivity, are calculated in the following way 

(NC)L. = (NC)S,, X ^ + (NC)£,, X (8) 

C' total total 

^ V ' " V ' 

(NC)|e (NC)fE 

where (NC)qe and (NC)^ are the NC backgrounds^ applicable for QE and IE events re- 
spectively. (Tqe, (Tie and Utotai are the relevant neutrino interaction cross-sections. 



C. Event Rates 

The number of (anti-)muon events^ in the i-th energy bin in the detector is given by 

^^ = irk / ■ dEA<P{E)a,^{E)RiE,EA)Pe,iE), (9) 



^ We assume that QE and IE events are fully uncorrelated with each other and that they can be cleanly 
separated. 

2 The (NC)^E and (NC)i^ backgrounds are smeared using the same energy resolution function that we use 
for QE and IE signal events respectively. Again, this choice is justified more by its apparent simplicity 
than by its realism. 

^ In principle, both detectors, also are sensitive to electron (positron) events. The number of electron 
events can be calculated using equation ^, by making appropriate changes to the oscillation probability 
and cross-sections. However, as was noted in [ p5[ , electron disappearance has hardly any sensitivity to 
013 and mass ordering at short baselines like the one under discussion. Also, it does not depend on Sep- 
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FIG. 2: Total event rates in five years as a function of sin^ 26'i3 in the FNAL - DUSEL set-up for -'^^Ne {ve) 
with 7 = 585 and ^He (t/g) with 7 = 350 are shown. The upper panels are for a 300 kt WC detector, while 
the lower ones are for a 50 kt LArTPC. Results are depicted for four different values of 5cp'- -90°, 0°, 90°, 
and 180°. Normal hierarchy has been assumed. For all other oscillation parameters we use the values given 
in equation nG. 



where T is the total running time, n„ is the number of target nucleons in the detector, 
e is the detector efficiency and R{E, Ea) is the GauBian energy resolution function of the 
detector. For muon (anti-muon) events, a^^ is the neutrino (anti-neutrino) interaction cross- 
section. The quantities E and Ea are the true and reconstructed (anti-)neutrino energies 
respectively and L is the baseline. 

We simulate the signal event spectrum using equation ^ for our assumed true"^ values for 



^ We distinguish between the "true" values of the oscillation parameters, which are used to compute the 
data, and their fitted values. Throughout this paper we denote the true value of a parameter by putting 
"(true)" after the symbol for the parameter. 



9 



the set of oscillation parameters as given in equation 0. The left-hand panels of figure |^ 
portray the number of neutrino events expected at DUSEL from five years exposure of the 
^'^Ne beta-beam from FNAL with 7 = 585 while the right-hand panels show the anti-neutrino 
events from five years run of the ^He beta-beam with 7 = 350. Results are presented as a 
function of sin^ 2^13 for four different values of 6cp- -90°, 0°, 90°, and 180° with normal mass 
hierarchy. The upper panels are for a 300 kt WC detector, while the lower panels are for a 
50 kt LArTPC. The number of events (both for ^^Ne and ^He ) varies in a wide range with 
the choice of 6cp- Out of the four different choices of 5cp, the maximum (minimum) number 
of events for neutrinos is obtained with 90° (-90°) irrespective of the choice of sin^ 26'i3. For 
anti- neutrinos, the same is true with 6cp — )■ Sep- We can explain this fact with the help of 
equation |l]. At the FNAL - DUSEL baseline, matter effects are small and hence A <^ 1 in 
equation |. The oscillation probabilities can be expressed as Pe^ — To+T_ sin 6cp+T^ cos 6cp 
and Pep, ~ To — T.sin^cp + T+cos6cp, where Tq,T± are independent of 6cp, whence the 
symmetry is manifest. For most of the energies T_ is positive assuming normal hierarchy. 
Now for 6cp = 90° (-90°) the term T_ sin(5cp gives a positive (negative) contribution towards 
the probability for neutrinos. The opposite is true for anti-neutrinos. The behavior of the 
number of events with sin^ 26'i3 is also understandable. As we increase the value of sin^ 2^13 
from a near- zero value, the 2nd and 3rd terms in equation |1|, which are linear in sin 2^13 and 
are dependent on 6cp, begin to contribute at first. Beyond a certain value of sin^2^i3, the 
first term, which does not depend on 6cp, takes over leading to the rise for all curves, but 
the relative size of CP effects decreases. This is the reason, why the discovery of leptonic 
CP violation does not become any easy task at large sin^ 2^13. 

As far as the total number of neutrino events is concerned, the performance of the different 
detectors is less different than indicated by a factor of six difference in the detector masses 
(see the upper and lower left panels of figure H). The reason is, that for the WC detector we 
consider only QE events but for the LArTPC we take into account both QE and IE events. 
At the same time the IE cross section is larger than the QE one at the relevant energies, see 
figure |. 



D. Reference set-up 



In order to compare the beta-beam set-up to alternative possibilities at FNAL, we in- 
troduce the wide-band beam concept. This project has been studied in detail in ^ . 
Here, a conventional neutrino beam will be sent from FNAL to a 300 kt WC in DUSEL. 
This beam has a spectrum wide enough to cover the first and second oscillation maximum 
and therefore can resolve most degeneracies. We use the same implementation as in |p3 |, 
with exception of the beam spectrum, which has been updated to correctly represent the 
beam which would be produced at FNAL |56|. Also, the beam intensity now corresponds 
to 1.2 MW, this is the maximum beam power which can be achieved at FNAL without 
Project X ||5^. The resulting physics sensitivities have been computed using GLoBES 



58 



and are shown as green shaded regions in figure I 
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5. THE NUMERICAL TECHNIQUE 



For all calculations we use as central (true) values 



Am2i 



2.4- 10-^ eV^ ±5%^ 
7.6-10"^eV2±2%, 



sin2 2e23 = 1.0 ±1% 
sin^ 9i2 = 0.32 ± 6% . 



(10) 



In all fits, these parameters are allowed to vary within the stated 1 a intervals. The central 



values are the current best fit values [59|; also the errors on the solar parameters are taken 
from The errors on the atmospheric parameters correspond to the result which are 

expected from T2K and NOz/A 

Here, we describe in detail the numerical procedure adopted to calculate the discovery 
potential of the FNAL - DUSEL beta-beam set-up. For our statistical analysis we use the 
following functions for WC detector and LArTPC 



total mc 



_ 2 _|_ 2 

~ X{ue^V^)Q-E ' X(i7e 

~^ X{Ue-^Ue)QE ~^ X(l7^ 

+ x" ■ 

' A,prior 



and 



{X 



totai J LArTPC 



The is given by 



i=l 



Xi 



where n is the total number of bins and 



(11) 



X{u^^U^)qE "t" X{lJe^l'f_i)lE X{ue^Ufj,)QE "I" X{Ue^l'^)i 
X(ue^Ue)QE ~^ X(ue^Ue)lE ~^ X{Ue^Ue)QE ~^ X(ue^l7e)lE 

+ x^ ■ . 

' A,prior 



(12) 



(13) 



U{^}. {es,a}) = Nr\{uj}) [1 + r^^Q + iVf 1 + A 



(14) 



Above, Nl^{{uj}) is the predicted number of QE events (calculated using equation |) in the 
z-th energy bin for a set of oscillation parameters u and are the number of background 
events in bin i. The quantities vr'^ and vr'' in equation |14| are the systematical errors on 
signals and backgrounds respectively. We consider tt^ = 2.5% and vr'' = 5% (see table 0). 
The quantities and are the pulls due to the systematical error on signal and background 
respectively. The data from equation O enters through the variable Xi = Nf^ + iVj*, where 



Nf^ is the number of observed QE signal events in the detector and is the background, 
as mentioned earlier. We simulate the QE signal event spectrum using equation ^ for our 



assumed true values for the set of oscillation parameters given in equation |T0|. We consider 
all the values of sin^ 26'i3 (true) and 5cp (true) in their allowed range and assume NH as 
true hierarchy. In a similar way, we estimate the contributions towards x^otai coming from 
other oscillation channels and event types (for both neutrino and anti-neutrino modes). In 



11 



CP violation discovery, NH true 





1 so 




IOC 

135 




90 




45 











Oh 




u 




CO 


-45 




-90 




-135 




-180 




10~ 



-l 1 — I — I I I I 



3a 



i i 



WC:A - 
WC:B - 



-3 



-2 



10 



10 



10 



sin 2013 ^^^^^^ 



FIG. 3: CP violation discovery potential of beta-beam at 730 km as range of Sep (true) as a function 
of the sin^ 2di3 (true) assuming normal hierarchy as true hierarchy at the 1 d.o.f. 3a C.L. Here we take 
7 = 350 for both ^^Ne & ^He ions. The results are shown for 300 kt WC detector assuming two different 
simulation methods. See the text for details. 



our fit we marginalize over all oscillation parameters and as well as the neutrino mass 
ordering, as applicable. We perform this by allowing all of these to vary freely in the fit 
and picking the smallest value of the function. However, we assume that some of these 
parameters which are poorly constrained by this experimental set-up, will be measured 
better from other experiments. Therefore, we impose a prior, or external constraint, on 
these parameters through Xprior^ given by 



|Am|;^| — (true) 



prior 



a(|Am2 



+ 



/ sin^ 26123 - sin^ 26*23 (true) ~ 



311 



+ 



' Am2i — Amgi (true) 
a(Amii) 



V a(sin" 2^23 

/ sin^ 6*12 — sin^ 9 
\ a(sin2 6*12 



'12 



(true) 



(15) 



where the la errors on these, are given in equation |T^. We minimize the Xtotai using the 



same procedure as it was described in the appendix of |23 . 



6. RESULTS & SUMMARY 



We evaluate the physics reach of the FNAL - DUSEL beta-beam set-up in terms of its 
discovery potentials for sin^ 26*13, CP violation and the mass hierarchy. These discovery po- 
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tentials quantify for any given sin^ 2^13 (true) for which range of possible values of Sep (true) 
the corresponding quantity will be discovered or measured at the chosen confidence level. 
The discovery reach for sin^ 2^13 is defined by the minimum value of sin^ 26*13 (true) which 
allows us to rule out sin^ 2^13 = in the fit. The CP violation discovery potential is defined 
as the range of 6cp (true) as a function of sin^ 2^13 (true) for which one can use the data to 
exclude the CP conserving solutions 6cp = 0° and 6cp = 180°. The mass hierarchy discovery 
reach is the limiting value of sin^ 2^13 (true) for which the wrong hierarchy can be excluded. 

Before we present our results for FNAL - DUSEL beta-beam set-up, we would like to 
discuss the CP violation discovery potential of beta-beam using a 300 kt WC detector at 
730 km taking 7 = 350 for both ^®Ne & ^He ions (see figure ^). This setup has already 



been considered in the literature and has been shown provide high performance The 



first and second oscillation maxima for 730 km baseline are at 1.4 GeV and 0.5 GeV for 
Am|^ = 2.4 ■ 10~^eV^. At this baseline, both the and the beam peak at the first 
oscillation maximum for 7 = 350. Also for this setup, we use two methods to treat the NC 
background in the water Cherenkov detector as described in section [4 A| . For method A, we 
have a threshold energy of 1 GeV for both ^^Ne and ®He ions and assume no backgrounds 
above the threshold. Since there is no background left above 1 GeV, it is possible to probe 
CP violation at very small values of sin^ 26*13 (true) as depicted by the dash-dotted blue 
line of figure ^ The major drawback of this method is that we do not have any access 
to the second oscillation maximum. Therefore, the {sgn{Am'^i),Scp) degeneracy for large 
values of sin^ 26*13 (true) is fully developed; as a result there is a sizable gap in sensitivity 
around sin^ 26*13 (true) = 3- 10~^. This effect has been described in \ WL\ and has been termed 
vr-transit. 

For method B, shown as dash-double-dotted black line in figure 0, where the threshold 
is 0.2 GeV and a background rejection factor of 10"'^ is used for both the ions. Although 
the total backgrounds are higher, the second oscillation maximum can be used, nonetheless. 
Therefore, the effects of vr-transit are mitigated for most of the parameter space, i.e. the 
sensitivity at large ^13 has essentially no gaps. The higher backgrounds result, however, in a 
reduced reach for small values of ^13. This concludes our comparison with previous results. 

Our results are summarized in figure ^, where the physics sensitivities for sin^ 26*13 (left 
panel), CP violation (middle panel), and for the mass hierarchy (right panel) are shown for a 
baseline of 1300 km. Results are shown at the 3 a confidence level for one degree of freedom. 
The various line styles are for the different detector options and background levels as given 
in the legend. The green shaded areas are the corresponding results for a conventional 



superbeam from FNAL to DUSEL using a 300 kt WC detector as described in detail in p5 



For the WC detector, we observe a distinct difference in sensitivities between the two dif- 
ferent schemes to include background (black, dash-double-dotted line and blue, dash-dotted 
line). Clearly, using a hard energy cut (WC: A, blue, dash-dotted line), which eliminates es- 
sentially all background and all information from the second oscillation maximum, performs 
better than WC: B (black, dash-double-dotted line). The only exception is found for the 
CP violation reach around sin^ 26*13 = 5 • 10~^ for a CP phase of approximately —120°. The 
reason, is the loss in information from the second oscillation maximum. Taking WC: A as 
our benchmark case, we find that a beta beam would result in about one order of magnitude 
improvement over a superbeam for all three measurements for at least one half of the CP 
phases. The performance of a 50 kt LArTPC is rather similar to the one of a six times larger 
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FIG. 4: Performance of FNAL - DUSEL beta-beam set-up at the 1 d.o.f. Scr C.L. in addressing 
sin^ 26113 , CP violation, and mass hierarchy discovery potential as range of ^cp (true) as a function of 
the sin^ 29i3 (true) assuming normal hierarchy as true hierarchy. The results are shown for 300 kt WC 
detector (assuming the simulation methods WC: A & WC: B. See Section 4 for details) and 50 kt LArTPC 
(assuming the background rejection factor of 10^^ and 10^*). The green shaded regions show the sensitivity 
of a wide band beam using the WC detector, as defined in detail in section 4D. 



WC detector using background scheme B. 

In comparison to the superbeam, shown as green shaded regions and described in sec- 
tion 4D, we see that using the same 300 kt WC detector the gain from using a beta-beam 
is obvious and most pronounced for < ^^p < 180°. This statement holds also for the use 
of a six times smaller 50 kt LArTPC. It is left to the judgement of the reader, whether the 
physics gain from a beta beam is commensurate to the effort of building a decay ring and 
of running the Tevatron for another decade. We also have shown that, given the perfor- 
mance parameters of the Tevatron, a baseline of around 730 km would provide a much more 
compelling physics gain. However, this baseline violates the boundary condition of locating 
the far detector at DUSEL. In summary, while a beta beam is a very interesting option to 
pursue high precision neutrino physics, especially the search for CP violation, it seems not 
to fit very well with existing and planed infrastructure in the US. 
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